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Abstract 

The KLOE experiment has collected 2.4 pb -1 of integrated luminosity during 
the commissioning of the DA$NE ^-factory in 1999. The performance of 
the detector has been studied using the <j> — » K^Ks decays collected during 
this period, yielding also first measurements of relevant K parameters such as 
masses and lifetimes. A clean K$ — * ir + ir~ sample is used to select Kl — * 
7T + 7r _ CP-violating decays and Kl — > Ks regeneration events in the detector 
material. Results on the regeneration probability in a beryllium-alluminum 
alloy and carbon-fiber plus aluminum composite are presented. 
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1 Introduction 



The DA$NE factory |J has come into operation in july 1999, delivering 2.4 
pb~ l of integrated luminosity to the KLOE experiment during its commissioning 
period. Such data allowed to perform a careful check of the detector operation and 
performances 0. 

A sample of — > K$Kl decays provided a source of monochromatic, 110 MeV/c 
{(3 = 0.21), beams of K$ and mesons, which were reconstructed in various decay 
modes. 

A preliminary measurement of the K$ lifetime was performed in a subsample 
of events. Charged decays oh the K L were also reconstructed, and a sample of 
CP violating Kl — > 7r + 7r~ decays was selected. The probability of Kl —> K$ 
regeneration in the materials of the KLOE detector was also measured. With the 
available statistics it was possible to clearly identify regeneration events produced in 
the beam-pipe, made of a Beryllium- Aluminum alloy, and in the drift chamber wall, 
made of Carbon fiber, and to measure the angular dependence of the cross section. 

2 Experimental setup 

The KLOE experiment |2| consists of a large volume tracking detector surrounded by 
a hermetic calorimeter both immersed in the axial magnetic field of 0.56 T produced 
by a superconducting coil. The tracking detector is a cylindrical drift chamber with 
alternated stereo views. The sensitive volume, 25 cm inner radius, 195 cm outer 
radius and 320 cm length is filled with a 90%Helium-10%Isobutane gas mixture. 
Charged particle tracks are reconstructed in 12582 drift cells arranged in 58 concen- 
tric layers. The spatial resolution in the transverse plane is about 150 jum. 

The calorimeter, made with 1 mm diameter scintillating fibers immersed in 0.5 
mm grooved lead plates, has a thickness of 23 cm corresponding to 15 Xq. It is 
segmented in 288 sectors in the barrel and 200 in the endcaps, each sector being 
subdivided in five longitudinal samplings. Energy clusters are reconstructed with 
a resolution cte/E = 0.057/ [E (GeV)] 1 ^ 2 © 0.006 and a space resolution of few cm. 
Time of flight measurement with a resolution of a t = 0.05 ns/[E (GeV)} 1 ^ 2 ©0.14 ns 
is used for particle identification. 

Data were taken at the energy of the resonance for six weeks in the period 
august-december 1999 during the commissiong of the DA$NE collider. Typical 
values of the beam parameters are listed in Table [|. The beams cross with a period 
multiple of 2.7 ns and with an angle of 25 mrad in the horizontal plane. This implies 
a small boost of 13 MeV of the center of mass relative to the experiment. 

numer of bunches 20 e + + 20 e~ 

beam currents 300 mA 

beam luminosity 10 30 cm~ 2 s~ l 
luminosity lifetime ~ 30 min 

Table 1: Typical values of the DA$NE beam parameters during data dating. 

The trigger requirement during most of the data taking was at least two calorime- 
ter energy deposits in the configurations Barrel-Barrel or Barrel- EndCap with E(Barrel)> 
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50 MeV and E(Endcap)> 90 MeV. This trigger has an average efficiency || of about 
90% for most of the decays and of about 84% for <fi — > K^Ks events with Kl and 
Ks decaying into charged particles of interest for this analysis. The trigger rate was 
typically 1.5 kHz mostly due to cosmic rays. The luminosity was measured with 
an accuracy of about 5% by recording Bhabha scattering events in the polar angle 
interval 22° < 9 < 158°. 

2.1 Data sample 

About 8 million <p decays were collected during data taking. The data were filtered 
against cosmics and machine background, and then classified in 5 different classes 

I: 

1. -> K S K L 

2. (j) -> K+K- 

3. 4> radiative decays 

4. pit events 

5. Bhabhas 

This analysis used only events in class 1 where a candidate Ks charged decay was 
present, satisfying the following requirements: one vertex made with two tracks of 
opposite curvature with p = (x 2 + y 2 ) 1 ^ 2 < 4 cm, \z\ < 8 cm, two pion invariant 
mass 400 < m < 600 MeV, and total momentum 50 < p < 120 MeV. 

3 Measurement of the Ks lifetime 
3.1 Measurement Technique 

The position of the secondary vertex (SV) of 7r + 7r~ pairs from K s decays, together 
with the knowledge of the position of e + e~ primary vertex (PV, the luminous re- 
gion) allows the measurement of the K s decay length in the frame, X s = (3^ct s . 
Xpy and Ypy are reconstructed run-by-run from Bhabha scattering events with a 
typical accuracy of few tens of microns and have widths L(X), L(Y) which are small 
compared to the 6 mm K s decay length ||J], ||. L(X) = 1.0 mm is measured by 
KLOE and L(Y) = few tens of microns is measured by DA$NE. Zpy is also re- 
constructed run-by-run with 100-200 fim accuracy, but it has a width of 12-14 mm, 
which distorts the K s lifetime distributions. Therefore, Zpy is computed event by 
event using the polar angle (6 S ) of the K s momentum (p s ) : 

Zpy = Z S y + COt(6 s ) X ^(X P y - X SV ) 2 + (Y P y - Y S y) 2 . (1) 

The effectiveness of this procedure is demonstrated by extrapolating the K s 
trajectory from the SV to the PV along p s /p s . Fig.|l] shows the difference between 
the intercept and the PV coordinate for x, y, and z, respectively (the selection of 
the data sample in these plots is described in the following section). 
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Figure 1: Projections of p s distance of closest approach to the PV 
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To avoid the systematic distortion of the lifetime distribution due to the SV 
resolution (which from Figs.[l] is shown to be comparable to the K s lifetime itself), 
the following estimator is used: 

A = if sv - f PV ) ■—. (2) 

Ps 

A is the K s decay length projected on the direction. A s can then be extracted 
from the fit to the A distribution using an exponential appropriately smeared with 
PV and SV resolution functions. Since from direct measurement L(X) = 1.0 mm 
and from Monte Carlo studies a(6 s ),a((p s ) ~ few tens of mrad, the resolution on 
the plots of Fig.Q can be interpreted in terms of the SV resolution, that turns out 
to be by far the dominant contribution. 

Finally, A is Lorentz-transformed to the <j) system using the run-average <fi boost 
as measured from Bhabha events |J. 

3.2 Sample Selection 

This analysis was performed using only a subsample of events for which the position 
of the luminous region had been computed after data taking. The selection required 
the simultaneous presence of the Ks candidate vertex and of a second two-track 
vertex outside a sphere of 11 cm radius. If there is a third vertex the event is 
rejected. The following cuts are then imposed to the K s candidates inside the beam 
pipe: 

(1) 493 MeV < M s < 497 MeV, 

(2) in the 0-CMS frame, 105 MeV < P s < 114 MeV, 

(3) in the 0-CMS frame, 45° < 6 S < 135°, 

(4) \cot9(ir,LAB)\ < 1.0 for both pions. 
The final selected sample is 6866 decays. 

3.3 Results 

The lifetime distribution for the 6866 decays was fit to four different functions, with 
the constraint of the total number of decays: 

(1) fit to an exponential smeared with two gaussian resolutions. 

(2) fit to an exponential smeared with three gaussian resolutions. 

(3) fit to an exponential smeared with two gaussian resolutions plus a third gaus- 
sian modeling a zero-lifetime component (due to background and/or to any 
other systematic effect). 

(4) fit to an exponential smeared with three gaussian resolutions plus a fourth 
gaussian modeling a zero-lifetime component (due to background and/or to 
any other systematic effect). 
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As an exanple cases (2) and (4) are shown in fig.0. In both fits the parameter PI 
is X s and P2 and P3 (P4 and P5) are the populations and standard deviation of 
the first (second) smearing gaussian. In case (2) P6 is the standard deviation of 
the third smearing gaussian, while in case (4) P6 and P7 are the population and 
standard deviation of the third smearing gaussian and P8 is the standard deviation 
of the fourth non-smearing gaussian. The results of the four fits are reported in 
table |. 



fit function 


X s (mm) 


2g smearing 

3g smearing 
2g smearing + lg non smearing 
3g smearing + lg non smearing 


5.71 ± 0.06 
5.71 ± 0.06 
5.90 ± 0.08 
5.83 ± 0.07 



Table 2: Typical values of the DA$NE beam parameters during data dating. 

The average X s value from the four fits is 5.78 mm. The fit statistical error on 
X s is taken as the largest fit error, 0.08 mm, while the fit systematic error is taken 
as half the difference between the maximum and minimum values of X s , 0.10 mm. 
The preliminary result is then: 

A s = 5.78 ± 0.08 (stat) ± 0.10 (syst) mm 

4 Measurement of the Kl — > K$ regeneration cross 
section 

Due to the different absorbtion cross section of K° and K°, a pure beam of Kl 
mesons will regenerate || Ks mesons when traversing the detector. Thus Kl — > 
Ks — * 7T7T is a potential source of background for the rare Kl — > 7T7t CP- violating 
decays in experiments where high precision in event counting is required. 

Regeneration is well studied at high energy but there is lack of experimental 
results at energies below 500 MeV. On the other hand theoretical predictions on the 
Kl — > Kg regeneration cross section at low energy suffer of large uncertainties. It 
has been shown [fjj that for a thin regenerator, t < Xs, the coherent regeneration 
is negligible. Inelastic regeneration can also be considered negligible in 110 MeV/c 
-fT^-nuclei interactions. Theoretical calculations || on incoherent Kl regeneration, 
based on the eikonal approximation and the Woods-Saxon form for the nuclear 
potential, foresee values in the range of 20-50 mb for the regeneration cross section 
on nuclei. 

4.1 Data analysis 

Ks — > 7r + 7r - candidates were selected as described in sect. 2.1. The additional 
requirement of the presence of at least one more vertex, with two unlike sign tracks, 
was made. 
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Ks — S» vtft lifetime fit: 3g smearing 




-4 -2 2 4 8 



DEMY LENGTH PROJECTED ON Ks MOMENTUM 



Ks ^ tt + 7t lifetime fit; 3g smearing, Ig prompt background 




DEMY LENGTH PROJECTED ON Ks MOMENTUM 



Figure 2: K s decay length fit to an exponential smeared with three gaussians without 
and with a fourth non-smearing gaussian. 
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Figure 3: Distribution of the Ks vertex invariant mass and momentum 



The distribution of the invariant mass and of the vector sum of the momenta in 
the (^-reference system are shown in FigfJ for the selected 604226 vertices (592687 
events). A clear peak at the Kg mass and at ps = 110 MeV/c is observed. The 
off-peak values are mainly due to n — > pi> decays in flight and to reconstruction 
errors. 

Fitting the distributions with two gaussians, the peak r.m.s. widths are a m = 1.0 
and a p = 1.9 MeV. 445347 well reconstructed K s — > n + 7i~ decays were selected 
requiring 

I m — (m) \ I p — (p) 



CT ... 




To search for the associated Kl decay vertex we first defined the Kl origin (the 
vertex) as the distance of closest approach between the Ks direction and the beam 
line (z axis). Events were retained requiring for the vertex position 

(£zMy + (»zMy + (£z«y< 9 

where a x = 0.63 cm, a y = 0.59 cm and a z = 1.38 cm are the r.m.s. widths of the <fi 
vertex distributions. 433786 vertices (433760 events) satisfied this cut. When more 
than one vertex were found in the Ks fiducial volume, the vertex associated with 
the invariant mass closer to the peak value was retained. This is the K$Kl sample. 

Associated charged decays of the Kl were then searched for. The K^ decay 
vertex is any second vertex reconstructed with two unlike sign tracks found in a cone 
opposite to the Ks direction in the reference system. Fig.[| shows the distribution 
of the angle 5 between the direction of the Kl, defined by —ps, and the line joining 
the Kl vertex and the <p vertex, for different intervals of the distance d between the 
two vertices. From such distributions we derived the angle resolution, as = 18 mrad, 
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Figure 4: Distribution of the angle between the Kl direction and the line joining 
the 0- vertex and the Kl- vertex for 5 < r < 15 cm, 25 < r < 35 cm, r > 45 cm. 



and the transverse vertex resolution, 5r± = 0.56 cm. The Kl vertex was accepted if 







crs 2 + 





1/2 



If more than one Kl vertex were found (only in 0.15% of the events) we selected 
the vertex with the smaller angle. This cut selected 134997 events. 

Different decays show different behaviour in terms of the missing momentum 
and the associated missing mass computed assuming the pion mass for all charged 
particles. Fig.^| shows the correlation of p mis and M^ is : K L — > 7r + 7r~7r° decays 
populate the region of M^ is = m^ ; K L — > n£v decays are clearly separated in two 
bands with M^ is < 0; K L — > 7r + 7r~ decays are peaked around p mis = 0, M^ is = 0; 
Kl — > Kg — > 7T + 7r _ "elastic" events are expected to populate the band with p mis = 
(_M 2 .) 1/2 - 

V mis/ 

Fig.[] shows the distribution of the Kl decay length. Fitting the distribution 
with an exponential in the region 40 cm < r < 150 cm, we obtain A = 333 ± 13 cm 
for the average Kl decay length in good agreement with the expected value of 343 
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cm, thus indicating a reconstruction efficiency constant over the range of the fit. 
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Figure 5: Correlation between the missing momentum and the squared missing mass 
for charged decays of the K L 



To select K L — > 7r + 7r~ decays and K L — > K s — > n + Ti~ elastic events we required 
the invariant mass associated with the Kl vertex to be the neutral kaon mass. Fig.[7] 
shows the invariant mass distribution fitted with a polynomial and a gaussian of 
r.m.s. width a m = 1.1 MeV. The distribution is divided in a signal band, m = 
497.7 ± 4 MeV, with 1991 events and two side-bands, m = 491.7 ± 2 MeV and 
m = 503. 7±2 MeV with a total of 1534 events populated mostly by K L semileptonic 
decays. 

The distribution of the decay distance is shown in Fig.^| for the signal band and 
the side-bands: the two peaks around r = 10 cm and r = 28 cm are interpreted 
as due to — > Ks regeneration in the spherical beam pipe and in the cylindrical 
drift chamber inner wall. 

For both Kl — > 7r + 7r~ decays and Kl — > Ks — > 7r + 7r~ elastic events the absolute 
values of the Ks and Kl momenta should be equal. Fig.|9| shows the distribution of 
the difference A[p| = \p$\ — \p^\ for the signal band and the side-bands. The signal 
band was fitted with a gaussian centred at A|p| =0 MeV with a = 2 MeV, repre- 
senting the decays and a second gaussian centred at A\p\ =3 MeV with o = 4 MeV, 
probably representing elastic interactions where a small fraction of the momentum is 
transferred to the recoil nucleus. To further reduce the background of semileptonic 
decays we required —6 < A\p\ < +12 MeV. This cut selected 915 events. 

Fig.0 shows the distribution of the vertex distance, r, and of the projected 
vertex distance, p = [x 2 + y 2 ) 1 ^ 2 , for the 915 events and for the side-bands sample. 
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Figure 6: Distribution of the K L decay length 
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Figure 7: Distribution of the invariant mass for charged decays of the Kl 
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Figure 8: Distribution of the decay distance for the signal and the side bands 




Figure 9: Distribution of the difference between the absolute values of the Ks and 
Kl momenta for the signal and the side bands 
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On the basis of the K$ decay length distribution of fig.@, two regions of interest were 
defined of width (-2, +4) cm around the position of the regenerators: 

• the beam pipe region, 8 cm < r < 14 cm, which contains 151 events 

• the chamber wall region, 23 cm < p < 29 cm, which contains 156 events. 

We parametrized the r and p distributions of the decay vertex using two gaussians 
for the peaks and a linear background, and we evaluated the amount of background 
events in the regions of interest. The number of Kl —>■ K$ — > 7r + 7r~ regenerated 
events is then: N% a = 123 ± 13, = 122 ± 12. 

The distribution of the angle, u, between the K$ and Kl momentum, shown in 



Fig.|TT|, was used to separate the K L — > decays, peaked at small angles, from 

the K L semileptonic decays and K L — > K s — > 7r + 7r~ elastic events. A cut at uj < 75 
mrad selected a sample of 279 Kl — > 7i + tt~ decays. 

4.2 Results 

The number of Kl — > Ks — > 7i + tt~ regenerated events is related to the regeneration 
cross section on nuclei, a reg , by 

N reg = N L E B S xn O reg 71 t (3) 

with Nl = Nsl e _r//A the number of Kl mesons reaching the regenerator, e the 
detection efficiency, Bsmr the Ks — > 7r + 7r _ branching fraction, n the number of 
nuclei of the regenerator per unit volume, t the thickness of the regenerator. 

The spherical beam pipe, of radius 10 cm, is made of AlBeMet, an alloy of 61%, 
in volume, of Beryllium (p = 1.85 g cm -3 ) and 39% of Aluminum (p = 2.7 g cm -3 ) 
and has a thickness of 0.50 mm: 

(nt) bp = 4.93 10 21 cm- 2 

The cylindrical drift chamber wall is made of Carbon 0.75 mm thick, 60%- fiber 
(p = 1.72 g cm -3 ) and 40%-epoxy (p = 1.25 g cm -3 ) and has a 0.20 mm thick 
Aluminum shield: 

(nt) cw = 6.97 10 21 cm- 2 

The crossing angle, computed event by event, gives an average increase of the 
thickness of 3% for the beam pipe and of 15.5% for the chamber wall. 

The detection efficiency for Kl — > Ks — > 7r + 7r~ regenerated events is the same 
as for Kl — > 7r + 7r~ decays selected in the analysis, the only difference being the final 
io cut. The efficiency is derived from the distribution of the Kl — * 7r + 7r _ events 
dNLmr/dr, the decay length in the laboratory A, and the number of KsKl events 

N SL B 

Ltttt —r fx I \ rfiV T l-k-k I . \ 

c e(r) = — (4) 

The experimental dNLnn/dr distribution for the 279 selected Kl — > 7r + 7r~ decays is 
shown in fig.|12|, where also the linear fit used to extract the results is shown. The 
region below 4 cm, coinciding with the Ks fiducial volume, has been excluded from 
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Figure 10: Distribution of the decay distance in space and of its projection on the 
x-y plane for K L — > -nil decays and regenerated events, and for the side bands sample 
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Figure 11: Distribution of the angle between the K$ and Kl momenta for the signal 
and the side bands 



the fit to avoid ambiguities, while the region above 140 cm has a drop in efficiency, 
being near the edge of the drift chamber volume. 
From the relations |3] and |] we obtain 

areg = X dN L /dr M W 

(dN ~L n7T i ' dr) was evaluated at the average radius of each region of interest, yield- 
ing: {dN^/drfP = 2.65 ± 0.19 cm" 1 , (dN L7Tn /dr) cw = 2.41 ± 0.15 cm" 1 . We 
obtain 

°%T M = ( 75 - 7 ± 9 - 6 stat) mb 

ag~ Al = (51.9 ± 6.2 stat ) mb 

where both the statistical errors coming from the event counting and the efficiency 
evaluation have been included. 

The angular distribution for regenerated events in the beam pipe and in the 



chamber wall is shown in Fig.|13|. To reduce the background only events in 10 cm 
< r < 13 cm, for the beam pipe, and in 24 cm < p < 27 cm, for the drift chamber 
wall were used. The small contamination from Kl — * tt + tt~ events is fully contained 
in the first bin, as was checked in the selected Kl — > 7r + 7r~ sample. The background 
from semileptonic decays is negligible, according to montecarlo estimations. 
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Figure 12: Radial distribution of the CP violating Kl — > vr + 7r decays, selected 
from data as explained in the text 



4.3 Systematic error 

Three main categories of systematic error sources have been considered: 

• event counting method; 

• efficiency evaluation; 

• regenerator thickness. 

For the first category two sources of error have been studied: the definition of the 
regions of interest (r.o.i) and the shape of the fit to the background. The r.o.i. limits 
were varied of ± 1 cm obtaining maximum variations of about 3% for the beam pipe 
events and of 2% for the drift chamber wall events. Various polynomial shapes were 
also fitted to the background, obtaining 2-=-3% variations on the beam pipe events, 
but much smaller (less than 1%) on the drift chamber wall. 

The efficiency evaluation is based on the assumption that regenerated K$ charged 
decays are detected with the same efficiency as the CP violating Kl — > tt + tt~ decays. 
A montecarlo study of such assumption [TI| shows that differences in efficiency can 
be expected up to 2%. 

A possible contamination of regenerated Ks decays in the Kl — > 7r + 7r~ sample 
could bias the measurement to higher efficiencies. To account for this effect the 



Kl — > 7r + 7r radial distribution, showed in fig. 12, was fitted adding to the main lin- 



ear shape two gaussians centered at the regenerators positions (just as was done in 
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fig.prOD . The data are compatible with absence of contamination on the drift chamber 
wall, but (probably due to a statistical fluctuation) allow a non negligible contami- 
nation on the beam pipe, which results in a considerably reduced efficiency (7-r8%). 
This turns out to be the main systematic error (non related to the regenerators 
knowledge) on this measurement. 

Finally the beam pipe thickness is known with a precision of 50 fim, accord- 
ing with the production tolerances. For the drift chamber wall an uncertainty of 
±50 fim on both carbon and aluminum was estimated from the study of the multiple 
scattering of e + e~ pairs. 

The systematic error sources are summarized in table [||. 



error source 


A(J Be-Al ( m6 ) 


Aa v ~ At (mb) 


r.o.i. limits 


2.5 


1.0 


backgr. shape 


2.0 


0.4 




1.5 


1.1 


regen. contam. in Kl — > 7t + tt~ 


6.5 


3.7 


regenerators thickness 


7.6 


3.5 


total 


10.6 


5.3 



Table 3: Systematic error contributions. 



5 Conclusions and discussion 

On the basis of the first data collected during the commissioning of DA$NE the 
KLOe experiment has analyzed ~ 1.4 10 5 — > KlKs decays with both kaons de- 
caying to charged particles. Using Bhabha scattering events to precisely define the 
collision region, and fitting the distribution of the K$ — > 7t + tt~ decay vertex, the 
Ks decay length is measured as: 

A s = 5.78 ± 0.08 (stat) ± 0.10 (syst) mm. 

The average Kl decay lenght is Xl = 333 ± 13 (stat) cm. A sample of 279 
Kl — ► 7T + 7r _ CP-violating decays is clearly identified with negligible background, 
and is used to measure the efficiency for reconstructing Kl — > Ks — > 7i + tt~ events 
due to regeneration in the beam pipe and in the drift chamber inner wall. 

The cross section in the beam pipe made of a 61%Beryllium-39%Aluminun alloy 

is 

ffBe-M = 75 7 ± g 6gtai ± 1Q mb 

The regeneration cross section in the drift chamber wall made mainly of Carbon is 

a°- Al = 51.9±6.2 stat ±5.3 syst mb 

With this data it is not possible to evaluate separately the regeneration cross section 
on Beryllium, Carbon and Aluminun nuclei unless we make additional hypotheses 
on its dependence upon the atomic mass. Theoretical calculations M based on the 



eikonal approximation are shown in Fig. 14 as well as a measurement on Beryllium 
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Figure 14: — » K$ regeneration cross section as a function of the atomic mass. 



ref.im H made with the CDM-2 detector at the Novosibirsk VEPP-2M electron- 
positron collider. The theoretical predictions do not show any evident dependence 
upon the atomic mass due to the interference of the K and K amplitudes that 
have different behaviour as a function of A. However our measurements can hardly 
accomodate a cross section on Alluminium nuclei smaller than for Beryllium and 
Carbon, as shown in fig.|15| where the correlation between the values of the cross 
sections is drawn. 
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